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Using Memory Tabu Search(MTS) algorithm, we investigate the relationship between structural 
characters and synchronizability of scale-free networks by maximizing and minimizing the ratio Q 
of the eigenvalues of the coupling matrix by edge-intercrossing procedures. The numerical results 
indicate that clustering coefficient C, maximal betweenness Bmax are two most important factors 
to scale-free network synchronizability, and assortative coefficient r and average distance D are the 
secondary ones. Moreover, the average degree (fc) affects the relationship between above structural 
characters and synchronizability of scale-free networks, and the minimal Q decreases when (fc) 
increases. 
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Based on nonlinear dynamics, synchronization in cou- 
pled dynamical systems has been studied for many years. 
It is observed in a variety of natural, social, physical, and 
biological systems and has found applications in a vari- 
ety of fields 0, S H S 0, H I3- In particular, synchro- 
nization in networks of coupled chaotic systems has re- 
ceive d a g reat deal of attention over the past two decades 
d, 0, irollTH. However, most of these works have been 
concentrated on networks with regular topological struc- 
tures such as chain s, g rids, lattices, and fully connected 
graphs [11 m El Il5l |. Recent empirical studies have 
demonstrated that many real-world networks have two 
common statistical characteristics: small- world effect [T^ 
and scale- free property 17], which cannot be treated as 
regular or random networks. Recently, an increasing 
number of studies have been devoted to investigating syn- 
chronization phenomena in complex networks with small- 
world and scale- free topologies l20> ■ 

One of the goals in studying network synchronization 
is to understand how the network topology affects the 
synchronizability. The network synchronizability can be 
measured well by the eigenratio Q of the largest eigen- 
value and the smallest nonzero eigenvalue p3 . l23 . l23.l25l| : 
thus, our work is to understand the relationship between 
network structure and its eigenvalues. Since there are 
several topological characters of scale- free networks, what 
is the most important factor by which the synchronizabil- 
ity of the system is mainly determined? 

In this brief report, we studied the relationship be- 
tween structural characters and synchronizability of 
scale-free networks. Some detailed comparisons among 
various networks have been done, indicating the network 
synchronizability will be stronger with smaller hetero- 
geneity, which can be measured by the variance of degree 
distribution or betweenness distribution , but 

the strict and clear conclusions have not been achieved 
because that previous studies are of both varying aver- 



age distances and degree variances. Another extensively 
studied one is average distance D. Some works indi- 
cated the average distance D is one of the key factors 
to network synchronizability |28| . Howe ver^ we have not 
achieved the consistent conclusion [T^ I21I |^. Some 
researchers considered that the randomicity is the more 
intrinsic factors leading to better synchronizability [2^ . 
which means that the intrinsic reason making small- world 
and scale-free networks having better synchronizability 
than regular ones is their random structures. Recently, 
several researches examine the effect of clustering coeffi- 
cient on the synchronization by us ing Kuramoto model 
[sof or master stability function [sil Is^ . Other re- 
searchers focus on the role played by maximal between- 
ness Bmax , they found the network synchronizability will 
be better with smaller Bmax 12(11331. Zhao et. al '271 en- 
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hanced the synchronizability by structural perturbations, 
they found that maximal betweenness plays a main role 
in network synchronization [3^ . Zhou et. al pS] studied 
the average distance D to synchronizability by crossed 
double cycle. However, a network contain several sta- 
tistical characteristics, such as degree distribution -P(fc), 
average distance D, clustering coefficient C, maximal be- 
tweenness Bmax and so on. In the previous works, if one 
wants to show clearly how a structural character affects 
the network synchronizability, such as average distance 
D, he would investigate the network synchronizability 
with different D while keeping other structural charac- 
ters constant approximately. However, this method ne- 
glect the influence made by the initial network struc- 
tural characters. In fact, the network functions, such 
as the synchronizability, are affected by these character- 
istics simultaneously. Therefore, we should investigate 
these structural factors holistically. In order to find the 
real factors affect network synchronization and eliminate 
the influence made by the initial networks, we maximize 
and minimize the eigenratio Q by MTS algorithm [ssf 
from the same initial network. The structural characters 
which change dramatically from maximal Q to minimal 
Q are the key factors influence network synchronizability, 
while the ones change little is not. 
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We investigate the synchronizability of a class of 
continuous-time dynamical networks with scale-free 
topology. Based on the synchronization criterion, we 
maximize and minimize the ratio Q of the eigenvalues 
of the coupling matrix by edge-intercrossing procedures, 
which provides a way for observing the correlation be- 
tween the synchronizability and those characteristics by 
keeping the degree distribution unchanged. 

We start by considering a network of N linearly cou- 
pled identical oscillators. The equation of motion reads 

N 

= F(a;*)+a^GyH(a;J), i = l,---,7V, (l) 

where x = F(x) govern the local dynamics of the vector 
field in each node, H(x) is a linear vectorial function, 
a is the coupling strength, and G is a coupling matrix. 

Stability of the synchronous state x^{t) = x^{t), i = 
1, • • ■ , N can be accounted for by diagonalizing the lin- 
ear stability equation, yielding N blocks of the form 

~ JF(xs)Ci — aXill{(i), where J is the Jacobian op- 
erator. Replacing aXi by v in the equation, the master 
stability function(MSF) 22] fully accounts for linear sta- 
bility of the synchronization manifold. For a large class 
of oscillatory systems, the MSF is negative in a finite 
parameter interval I^t = (i^i < < 1^2) ^3- When the 
whole set of eigenvalues (multiplied by a) enters the inter- 
val 1st , the stability condition is satisfied. This is accom- 
plished when (TA2 > V2 and ctAat < 1^2 simultaneously. 
As 1^2 and depend on the specific choice of F(x) and 
H(x), the key quantity for assessing the synchronization 
of a network is the eigenratio 

Q = Ajv/A2, (2) 

which only depends on the topology of the network. The 
small Ajv/A2 is, the more packed the eigenvalues of G 
are, leading to an enhanced a interval for which stabil- 
ity is obtained|36|. In this paper, we will not address a 
particular dynamical system, but concentrate on how the 
network topology affects eigenratio Q. 

The processes of heuristic algorithm, named MTS, is 
as follows. 

Step 1. Generate an initial matrix Go of the extensional 
BA network with N nodes and E edges. Set 

the optimal network' coupling matrix = Gq and 
the optimal network of taibu table Gk — Gq, and 
the time step fc = 0. Compute the ratio Q of G^. 

Step 2. If a prescribed termination condition is satisfied, 
stop; Otherwise intercrossing a pairs of edges cho- 
sen randomly based on the network remains con- 
nected, denote by G. 

Step 3. If the ratio Q of G, denoted by Qc, satisfying 
Qg < Qgi , Qgi = Qg, else if Qg < Qg^ , Gk+i := 
G. When Qg > QGk^ if does not satisfy the 
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FIG. 1: (Color online) The structural characters vs the eigen- 
ratio Q. (a) Average distance D. (b) Clustering coefRcient 
C. (c) Assortative coefficient r. (d) Maximal node between- 
ness Bmax- The blue dot line denotes the state of the initial 
network. The data are averaged over ten independent runs of 
network size A'^ = 500. 



tabu condition, IQg^ ~ Qg\/Rg > ^ (where (5 is a 
random number between 0.5 and 0.75), Gk+i = Gk, 
else Gfc+i = G. Go to Step 2. 

Since the MTS algorithm is heuristic, it can only find 
the approximate optimal solution. Thus, the termination 
condition of Step 2 should confirm by the experimenta- 
tion solution. 

The numerical results are experimented on extensional 
BA model for different network scales. The statistical 
properties of the optimal networks show similar trends. 
After many numerical experimentations, we set the ter- 
mination condition for maximizing Q as 8000 time steps 
and the one for minimizing Q as 3000 time steps, which 
can obtain the stability value using MTS algorithm. 

We start from a network of size N = 100, 200, 300, 
400, 500 and the average degree (fc) = 6 and then per- 
form the optimization precesses. At each time step, we 
record the structural properties, such as D, G, r and 
average node betweenness, when the objective function 
Q is reduced. Let D™", G"™ r"™ and S™™ denote 
the stable value when Q reaches its minimal value Qmim 
and i:)™'^^, G™°^, r"'^^ and B^^i^^ denote the stable value 
when Q reaches its maximal value Qmax- Define the rel- 
ative diversity function of the structural character x as 
follows 



which can denote the difference of structural character 
X to Qmax and Qmin- The larger f{x) is, the structural 
character x change dramatically when the network leave 
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FIG. 2: (Color online) Bmax obtained by maximizing and 
minimizing Q vs network size A*'. 
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FIG. 3: (Color online) Relationship between Q and average 
degree (fc). The inset shows the functions of the structural 
characters to different (fc). 

far from its optimal synchronizability state, which means 
X is more relevant to network synchronizability. 

Figure^ (a) demonstrates that D remains stable when 
minimizing Q, while increase a little when maximizing Q. 



Figure^ (b), (c) show C and r decrease to a stable value 
when maximizing and minimizing Q, and C"""^ and r™*^^ 
are both smaller than (7™™ and r"*™ . The difference be- 
tween C™" and (7™"^, r™™ and r™"^ means that the 
two structural characters are relevant to synchronizabil- 
ity of scale- free networks. Figure ^ (d) gives the change 
trend of average node betweenness when maximizing and 
minimizing Q, which is consistent with D. Figure |21 (a)- 
(d) demonstrate the stable value of _D, C, r and Bmax 
when N = 100,200,300,400,500. From Fig|21 one can 
obtain that when N ^ 500, f{D) = 1.27, /(C) = 20.97 
and /(r) = 1.27 and f{B„iax) = 19.04. Moreover, one 
can see that the relative diversity of C, r and Bmax be- 
come large, while the one of D remain constant, which 
indicates that the influence produced by the structural 
characters C, r and Bmax to synchronizability of scale- 
free networks would become great when TV become large. 
Furthermore, we investigate the relationship between av- 
erage degree (fc) and f{x). Figure O dcmonstrats the 
Qmin obtained by MTS algorithm to different (k) when 
N = 500. The inset gives the functions of the structural 
characters obtained by different (k). From the inset, one 
can see that if {k) increase, the function of D, r and 
Bmax increases while the function of C decreases, which 
means that the influence of the structural characters to 
network synchronizability is affected by (k). When (k) 
increases, Bmax, D and r become more relevant to syn- 
chronizability of scale- free networks, while C become less 
relevant. 

In summary, using the MTS optimal algorithm, we 
maximized and minimized the network synchronizabil- 
ity by changing the connection pattern between differ- 
ent pairs of nodes while keeping the degree variance un- 
changed. Starting from cxtcnsional BA networks, we 
found the relationship between structural characters and 
synchronizability of scale-free networks. The numerical 
results indicate that D, C, r and Bmax influence net- 
work synchronizability simultaneously. Especially, C and 
Bmax are the two most important structural characters 
which affect synchronizability of scale-free networks, as- 
sortative coefficient r is the secondary character and D is 
the last one. Furthermore, the relationship is affected by 
the average degree (k) , and the maximal synchronizabil- 
ity of scale-free networks increases when (k) increases. 
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